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1.0  ABSTRACT 


The  wideband  coverage  of  the  traditional  Fast  Fourier  transfonn  (FFT)-based  electronic  warfare 
(EW)  receiver  is  limited  by  the  sampling  frequency  of  the  analog-to-digital  converter  (ADC). 
Compressed  sensing  (CS)  based  wideband  EW  digital  receiver  has  been  proposed  and  under 
development  to  process  wideband  sparse  signals  with  reduced  sampling  rate.  In  this  study,  we 
calculate  the  sensitivity  of  a  CS  based  EW  receiver.  Two  modulation  schemes  of  CS  are  used. 
One  is  pseudorandom  chip  (PRC)  and  the  other  is  non-uniform  sampling  (NUS)  modulation. 

The  orthogonal  matching  pursuit  (OMP)  algorithm  is  used  to  process  CS  data.  The  insufficient 
nature  of  the  sparcity  of  the  signal  adversely  affects  the  signal  detection  probability  for  both  PRC 
and  NUS  modulation  schemes.  Methods  to  remedy  this  problem  in  both  modulation  schemes  are 
proposed.  The  improvement  in  the  signal  detection  and  thus  the  sensitivity  perfonnance  of  the 
CS  based  receiver  are  presented. 
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2.0  INTRODUCTION 


Electronic  Intelligence  Receiver  (ELINT)  is  an  important  component  in  EW  and  layer  sensing. 
The  infonnation  provided  by  its  function  of  constant  surveillance  can  be  used  to  detect,  track  and 
classify  signals  across  the  electromagnetic  spectrum.  The  proper  identification  and  reaction  to 
the  threat  can  avoid  disaster  and  assure  spectrum  dominance  for  Air  Force  systems.  To  meet  the 
challenges  in  today’s  and  tomorrow’s  EW  environment  where  signals  are  increasingly 
sophisticated,  more  dynamic  and  more  crowded  in  radio  frequency  (RF)  spectrum  compared  to 
yesterday,  digital  wideband  receiver  is  under  intensive  development.  The  current  state-of  the-art 
for  digital  receiver  bandwidth  coverage  is  now  reaching  multi-GHz.  The  conventional  wideband 
digital  ELINT  receiver  design  is  based  on  the  Nyquist  infonnation  theory,  whereby,  bandwidth 
coverage  is  limited  by  the  Nyquist  sampling  rate.  Perfonnance  therefore,  depends  on  the  current 
high  speed  ADC  technology  and  computation  hardware  such  as  field-programmable  gate  array 
(FPGA).  While  these  technologies  are  advancing,  the  progress  of  the  improvement  may  not  keep 
up  with  the  demand  of  ELINT  requirements  driven  by  next  generation  cognitive  and  broadband 
EW  capability  needs.  Innovative  designs  are  proposed  to  keep  abreast  with  the  requirement 
while  to  stay  within  the  current  technology  limitation.  One  potential  leap-ahead  is  to  apply 
compressed  sensing  concept  for  ELINT  receiver  development. 

CS  technique  provides  an  alternate  solution.  It  has  received  significant  research  attention  since 
the  fundamental  theory  was  proved  by  E.  Candes,  J.  Romberg,  and  T.  Tao  in  2006  [1].  In  this 
paper,  the  authors  indicated  that  the  band  limited  by  the  Nyquist  sampling  rate  is  too  restrictive. 

It  is  particularly  true  when  the  signal  is  sparse,  i.e.,  the  signal  can  be  represented  by  only  a  few 
significant  components  in  an  orthogonal  basis.  Under  the  sparse  condition,  one  can  devise  a 
sampling  scheme  that  takes  a  much  reduced  data  rate  compared  to  that  of  Nyquist  sampling  rate 
and  still  cover  the  same  bandwidth  as  the  Nyquist  sampling  case  does.  The  sampling  scheme 
usually  involves  modulating  the  incoming  signal  and  sampling  a  data  after  summing  the 
modulated  signal  over  an  interval  of  several  chips.  The  Nyquist  waveform  can  be  reconstructed 
using  the  reduced  sampling  if  the  measurement  matrix  representing  the  sampling  scheme  has  the 
restrictive  isometric  property  (RIP). 

In  this  study,  we  include  two  measurement  matrices  which  both  exhibit  high  probability  of  RIP. 
One  is  the  PRC,  and  the  other  is  NUS.  The  PRC’s  measurement  matrix  element  is  a 
pseudorandom  bit  of  +  1  and  -1  with  equal  probability.  The  scheme  was  proposed  [2]  for  image 
processing  using  single  pixel  camera,  where  the  field  of  view  was  masked  by  a  grid  structure  of 
random  +1  and  -1  pattern.  The  NUS  scheme  was  proposed  and  fabricated  by  Defense  Advanced 
Research  Projects  Agency’s  (DARPA)  analog-to-information  (A-to-I)  program  initiated  in  2008 
[3].  In  this  scheme,  an  internal  Nyquist  rate  (i.e.,  chip  rate  of  modulation)  of  4.8  GHz  is  applied. 
The  sampling  time  is  a  pseudorandom  integer  number  within  a  certain  range  of  time  grid 
determined  by  the  internal  Nyquist  rate.  The  average  sampling  rate  is  about  300  MHz  and  the 
average  down  sampling  rate  is  about  16. 

Many  reconstruction  algorithms  have  been  proposed  to  restore  the  original  Nyquist  wavefonn  of 
sparse  signal  from  the  reduced  sampling  data  set.  These  reconstruction  algorithms  have  a  range 
of  complexity  and  accuracy.  We  apply  OMP  to  assess  the  sensitivity  of  the  receiver.  OMP  is  a 
restoration  algorithm  developed  at  an  early  stage  [4].  The  algorithm  is  to  locate  and  calculate  the 
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signal’s  component  iteratively.  There  are  two  major  computation  steps  in  every  iteration.  One 
is  to  correlate  the  residual  measurement  data  with  the  measurement  basis  to  locate  and  update  the 
signal  component.  The  other  is  to  update  signal’s  components  using  least  square  fit  scheme. 

The  computation,  compared  to  the  later  developed  algorithms  including  NESTA  [5],  is  less 
complex  and  yet  it  has  yielded  reasonably  accurate  frequency  results. 

2.1  CS-Based  Receiver  Sensitivity  Simulation  Using  OMP 

2.1.1  Routine  of  Receiver  Sensitivity  Simulation 

Receiver’s  sensitivity  calculation  deals  only  one  signal,  while  dynamic  range  (DR)  deals  two  or 
more  signals  [6],  The  DR  simulation  of  CS-based  receiver  will  be  conducted  in  future  study. 

There  are  two  stages  in  the  routine  of  sensitivity  simulation.  The  first  is  to  determine  the 
threshold  of  the  detection  which  depends  on  the  false  alarm  rate  [6].  In  this  stage,  one  needs  to 
build  an  ensemble  to  set  up  probability  density  function  (PDF)  of  the  stochastic  parameter  of  the 
observed  parameter  at  the  output  port  when  a  Gaussian  white  noise  is  at  the  input  port.  Result 
from  this  stage  is  the  detection  threshold  to  be  used  in  the  second  stage.  The  second  stage 
calculates  the  detection  probability  (Pd)  vs.  signal-to-noise  (S/N)  ratio.  The  default  false  alarm 
rate  is  1CT7  in  this  study.  At  each  S/N  point,  Pd  is  calculated  from  the  average  of  10,000  runs. 
Sensitivity  is  defined  at  the  S/N  (in  dB)  when  probability  of  detection  is  at  90%. 

2.1.2  Pseudo-Random-Chip  and  Non-Uniform  Measurement  Matrix: 

A  schematic  CS  representation  is  given  by: 


-v.t/xi  -  ^fi/x.vTVxi  -  ^i/x.y^.Vx.yEYxI  q 

where  y  and  x  are  the  measurement  and  Nyquist  vectors,  respectively.  The  ®  is  the 
measurement  matrix  with  a  dimension  of  M  by  N.  The  'P  is  the  basis  matrix.  The  Fourier  basis 
is  used  in  this  study.  The  s  is  a  sparse  vector  in  frequency  domain. 


A  schematic  of  ®  using  PRC  modulation  is  shown  as: 
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where  each  entry  in  the  matrix  is  either  +1  or  -1  with  an  equal  probability  among  them.  In  this 
study,  the  dimension  of  the  measurement  matrix  is  36x256,  and  thus  the  sampling  reduction  rate 
is  7.1. 
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A  schematic  of  ®  using  NUS  modulation  is  shown  as: 
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(3) 


where  each  row  has  only  one  non-zero  entry  of  1  and  the  rest  elements  are  zero.  The  distance 
between  the  adjacent  row’s  non-zero  entries  is  a  random  integer  number  within  a  rage.  In  this 
study,  it  is  between  7  and  26  with  a  unifonn  distribution.  The  dimension  of  the  NUS 
measurement  matrix  is  36x645  in  this  study,  and  thus  the  sampling  reduction  rate  is  17.9. 

2.1.3  Orthogonal  Matching  Pursuit 

The  OMP  algorithm  is  depicted  in  Figure  1  [5]: 


1: 

r<°> «-  y 

2: 

A<°>  «-  0 

3: 

for  i  —  1, . . . ,  k  do 

4: 

A(0  —  arg  max/=u  „  |(V'- 

5: 

aW  —  Ah"1)  U  Ah) 

6: 

q>(»)  <_  [oh- 1)  *AW] 

7: 

xh)  <—  arg  min*  ||y  —  Oh)x 

2 

8: 

ah)  <-  <J>h)xh) 

9: 

'’rl 

1 

1 

10: 

end  for 

11: 

x  <—  Xh9 

12: 

return  x,  A(<r),  aik\  r':k} 

Figure  1:  OMP  Algorithm 


where  steps  4  and  7  are  most  computation  heavy.  Step  4  is  to  find  signal’s  component  through 
correlating  the  residuals  of  the  measurement  vector  with  each  column  of  the  measurement 
matrix.  Step  7  is  to  calculate  the  updated  Nyquist  vector  through  least  square  method.  Since  the 
incoming  signal  is  presented  using  real  number  and  the  Fourier  basis  matrix  is  a  complex  matrix, 
there  are  two  components  in  Fourier  bases  for  each  frequency  components  in  the  incoming 
signal. 
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2.2  Results  of  CS-Based  Receiver  Sensitivity  Simulation  Using  OMP 
2.2.1  PRC  Results  and  Discussion 

Figure  2  shows  the  PDF  of  the  PRC  CS-based  Receiver  using  OMP.  The  PDF  can  be 
approximated  with  a  two  half  Gaussians  and  a  fat  tail  (i.e.,  exponential  decade)  as  in  the 
following  equation: 

-  (x  -  0.3161)2 

p(x)  =  16898exp[ — - — ];  0.3161  >  x 

0.004706 

,  -  (x  -  0.3161)2 

3  p(x)  =  16898exp[ — - —];  0.5233  >  x  >  0.3161  (4) 

0.03302 

p{x)  =  4636exp[-20.82(x  -  0.5233)];  x  >  0.5233 

where  x  is  the  random  number  of  the  output  of  the  PRC-based  receiver  using  OMP  when 
Gaussian  noise  is  at  the  input.  The  threshold  is  1.159  when  false  alarm  rate  is  10’7. 


Figure  2:  Probability  Density  Function  of  PRC  CS  Scheme  Receiver 

Figure  3  shows  the  detection  probability  vs.  S/N.  There  are  two  curves  in  this  figure.  One  is  the 
result  according  to  the  detection  threshold,  and  the  other  is  the  detection  according  to  the 
accuracy  of  the  frequency  determination.  For  the  latter,  the  detection  is  claimed  if  the 
determined  frequency  is  accurate  within  one  Fourier  frequency  bin  (frequency  resolution  in  FFT 
of  the  frame-sized  data  points).  One  notes  that  the  frequency-determined  detection  only  reach  to 
about  96%  when  S/N  increases.  Its  significance  is  that  there  is  about  4%  that  frequency  is 
wrongly  determined  even  signal  is  claimed  according  to  threshold  criterion.  This  nuisance 
deserves  attention.  A  close  look  into  cases  giving  wrong  frequency  results  found  the  problem  is 
due  to  sparcity.  This  is  evidenced  when  the  input  frequency  is  not  at  the  Fourier  bins,  and  the 
signal  power  spreads  over  to  the  neighboring  bins.  The  default  rectangular  window  results  in  a  - 
13  dB  sidelobe,  which  may  be  too  big  for  OMP  to  yield  accurate  solutions. 
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Figure  3:  Detection  Probability  vs.  S/N  for  PRC  Case  without  Remedy 

The  measurement  matrix  size  is  36x256.  *:  detection  based  on  detection  threshold;  o:  detection 

based  on  frequency  accuracy. 


A  simple  remedy  for  this  problem  is  to  apply  Blackman  window  which  is  known  to  reduce 
sidelobe  to  -58  dB.  Repeating  the  sensitivity  calculation,  the  result  after  this  remedy  is  shown  in 
Figure  4.  Here  only  the  Pd  vs.  S/N  is  shown  since  the  PDF  result  is  similar  to  Figure  2  and  not 
shown  here.  One  notes  that  the  gap  disappears  and  the  detection  probability  reach  nearly  100% 
as  S/N  increases.  One  also  notes  that  the  sensitivity  is  about  6.8  dB  according  to  threshold 
detection  (1  dB  according  to  frequency-determined  criterion). 


Figure  4:  Detection  Probability  vs.  S/N  for  PRC  Case  with  Remedy 

The  measurement  matrix  size  is  36x256.  *:  detection  based  on  detection  threshold;  o:  detection 

based  on  frequency  accuracy. 
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2.2.2  NUS  Results  and  Discussion 


Figure  5  shows  the  PDF  of  the  NUS  CS-based  receiver  using  OMP.  Same  as  the  PRC  case,  the 
PDF  can  be  approximated  with  a  two  half  Gaussians  and  a  fat  tail  (i.e.,  exponential  decade)  as  in 
the  following  equation: 

—  (jc  —  0.3834)2 

p(x)  =  19141exp[ — - — ];  0.3834  >  x 

0.003930 

-  (x  -  0.3834)2 

p(x)  =  19141exp[ — - — ];  0.6168  >  x  >  0.3834  (5) 

0.012776 

p(x)  =  272exp[-26.99(x  -  0.6168)];  x  >  0.6168 


where  x  is  the  random  number  of  the  output  of  the  NUS-based  receiver  using  OMP  when 
Gaussian  noise  is  at  the  input.  The  threshold  is  1.003  when  false  alarm  rate  is  10'7.  Figure  6  is 
the  detection  probability  vs.  S/N. 


Figure  5:  Probability  Density  Function  of  NUS  CS  Scheme  Receiver 


Figure  6:  Detection  Probability  vs.  S/N  for  NUS  Case  without  Remedy 

The  measurement  matrix  size  is  36x645.  *:  detection  based  on  detection  threshold;  o:  detection 

based  on  frequency  accuracy. 
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As  in  Figures  3  and  4,  the  two  curves  in  this  figure  are  the  detection  probability  using  either  the 
threshold  or  the  frequency  criterion.  Similar  to  Figure  3,  the  frequency-determined  detection 
does  not  reach  100%  when  S/N  increases.  First  thought  is  to  apply  Blackman  window  as  the  one 
used  for  PRC  case  to  remedy  the  problem.  However,  it  does  not  yield  improvement.  The  reason 
is  due  to  the  difference  in  NUS’s  measurement  matrix  which  is  not  as  robust  as  PRC’s.  Since 
there  is  only  one  non-zero  entry  for  each  row  in  NUS  measurement  matrix,  the  Blackman 
window  only  applies  to  that  single  non-zero  entry,  and  thus  exerts  little  effect  to  suppress 
sparcity.  A  close  look  into  those  cases  giving  wrong  frequency  results  are  shown  in  Figure  7.  It 
shows  that  when  input  frequency  is  at  the  Fourier  bins,  the  frequency  detection  reaches  nearly 
100%  as  S/N  increases.  However,  when  input  frequency  is  at  the  middle  of  two  adjacent  Fourier 
bins,  the  frequency  detection  reaches  only  80%  as  S/N  increases.  With  this  information,  we 
devise  a  remedy  for  NUS  CS  signal  detection. 


Figure  7:  Probability  of  Detection  vs.  S/N  when  the  Input  Frequency  is  restricted  to  at 
Quarter,  Half,  and  Zero  Integer  away  from  the  Fourier  Bins 

The  last  figure  has  the  input  frequency  randomized. 
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In  addition  to  the  computation  using  the  normal  Fourier  basis  at  the  integer  bins,  we  apply  an 
extra  computation  using  Fourier  basis  at  half  integer  bins.  The  matrix  and  matrix  element  for  the 
Fourier  basis  at  the  integer  bins  and  at  the  half  integer  bins  are  given  by: 


Fkn  =  exp  (-J2nkn/N  );  T  =  F 

F\/2,kn  =  exp  +  1 7  2)«/7V );  'T'1/2  =  FU2 


(6) 


where  'F  and  4^/2  are  the  matrices  for  Fourier  basis  at  integer  and  half  integer  bins,  respectively. 
By  comparing  the  results  using  T*  and  ¥1/2,  a  more  accurate  frequency  result  can  be  obtained. 


Repeating  the  sensitivity  calculation,  the  result  after  this  remedy  is  shown  in  Figure  8.  Here  only 
the  Pd  vs.  S/N  is  shown  since  the  PDF  result  is  similar  to  Figure  5  and  not  shown  here.  It  is 
noted  that  the  gap  disappears  and  the  detection  probability  reaches  nearly  100%  as  S/N  increases. 
One  also  notes  that  the  sensitivity  is  about  4.8  dB  according  to  threshold  detection  (0.8  dB 
according  to  frequency-determined  criterion). 


Figure  8:  Detection  Probability  vs.  S/N  for  NUS  Case  with  Remedy 

The  measurement  matrix  size  is  36x645.  *:  detection  based  on  detection  threshold;  o:  detection 

based  on  frequency  accuracy. 
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3.0  SUMMARY  AND  FUTURE  WORK 

Compressed  sensing  technology  has  generated  interest  in  image  and  signal  processing.  It  is  of 
great  interest  to  assess  the  perfonnance  of  CS-based  EW  receiver  in  terms  of  sensitivity  and 
dynamic  range.  We  have  conducted  sensitivity  simulation  of  the  CS  cases  using  PRC  and  NUS 
modulations  and  OMP  restoration  algorithm.  It  has  been  found  that  sparcity  poises  problem  in 
frequency  determination.  As  a  result,  the  frequency  cannot  be  accurately  calculated  even  the 
signal  is  detected.  For  PRC  case,  we  propose  applying  Blackman  window  in  the  measurement 
matrix  to  suppress  sidelobe  and  subsequently  to  improve  sensitivity  performance  of  the  PRC 
modulation  case.  For  NUS  case,  we  propose  an  algorithm  that  combines  Fourier  basis  at  the 
integer  bins  and  at  the  half  integer  bins  to  improve  the  frequency  accuracy.  The  results  of  the 
improvement  using  these  remedies  have  been  shown  in  this  study. 

The  future  work  will  continue  assessing  perfonnance  of  the  CS-based  receiver  using  different 
restoration  algorithms  including  NESTA.  It  will  also  include  the  hardware  implementation  of 
both  measurement  matrix  in  the  front-end  and  signal  processing  in  the  back-end.  The  front-end 
module  involves  circuit  design  to  realize  PRC  and  NUS  modulation  schemes.  The  back-end 
module  involves  loading  algorithm  into  computation  devices  such  as  general  purpose  computing 
on  graphics  processing  unit  (GPGPU)  for  signal  detection  and  classification  in  real  time. 
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LIST  OF  SYMBOLS,  ABBREVIATIONS,  AND  ACRONYMS 


ACRONYM 

DESCRIPTION 

ADC 

analog-to-digital  converter 

A-to-I 

analog-to-infonnation 

CS 

compressed  sensing 

DARPA 

Defense  Advanced  Research  Projects  Agency 

DR 

dynamic  range 

ELINT 

Electronic  Intelligence  Receiver 

EW 

electronic  warfare 

EFT 

Fast  Fourier  transform 

FPGA 

field-programmable  gate  array 

GPGPU 

general  purpose  computing  on  graphics  processing  unit 

NUS 

non-uniform  sampling 

OMP 

orthogonal  matching  pursuit 

Pd 

detection  probability 

PDF 

probability  density  function 

PRC 

pseudorandom  chip 

RF 

radio  frequency 

RIP 

restrictive  isometric  property 

S/N 

signal-to-noise 
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